# To excute:
# ?> cat readnplot.r | rloc --slave &
#
#load(".RData")
# Load functions
source("../readf.r")
source("../plotfunc.r")
source("../sigtest.r")

# Read in experiment output
print("Read in output....")
#n10.0301 <- read0301("par.10.0301.p2",nexp=10)
#n10.0301 <- readf("par.10.0301",nexp=30)
#n10.0401 <- readf("par.10.0401",nexp=30)
#n10.0402 <- readf("par.10.0402",nexp=30)
#n10.0403 <- readf("par.10.0403",nexp=30)
#n10.0405 <- readf("par.10.0405",nexp=30)
#n10.0406 <- readf("par.10.0406",nexp=30)
n10.0515 <- readf("par.10.0515",nexp=30)
n10.0516 <- readf("par.10.0516",nexp=30)
n10.0525 <- readf("par.10.0525",nexp=30)
n10.0526 <- readf("par.10.0526",nexp=30)
# organize data
print("Organize data")
#rm(pn10.objf,pn10.ofc)
rm(pn10a.objf,pn10a.ofc)
#pn10.objf <- list(at.as=n10.0301$objf, co.as=n10.0401$objf, co.ws=n10.0402$objf, co.ai=n10.0405$objf, co.wi=n10.0406$objf, co.mo=n10.0403$objf)
#pn10.ofc <- list(at.as=n10.0301$ofc, co.as=n10.0401$ofc, co.ws=n10.0402$ofc, co.ai=n10.0405$ofc, co.wi=n10.0406$ofc, co.mo=n10.0403$ofc)
pn10a.objf <- list(as.ai=n10.0515$objf, as.wi=n10.0516$objf, ws.ai=n10.0525$objf, ws.wi=n10.0526$objf)
pn10a.ofc <- list(as.ai=n10.0515$ofc, as.wi=n10.0516$ofc, ws.ai=n10.0525$ofc, ws.wi=n10.0526$ofc)
#print("")
# Plotting
#ploti(pn10.objf, natp=1024, nptp=256, ps=T, yrange=c(0.45,0.55))
#plotco(pn10.ofc, natp=1024, nptp=256, ps=T, yrange=c(-1,1))
# Clean up
#rm(n10.0401,n10.0402,n10.0403,n10.0405,n10.0406)
# 
#quit(save="yes")
